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Abstract. We show that the thermodynamics of the focusing cubic 
discrete nonlinear Schrodinger equation are exactly solvable in dimen- 
sions three and higher. A number of explicit formulas are derived. 



1. Introduction 

A complex- valued function u of two variables x and t, where x 6 M d 
is the space variable and t G M is the time variable, is said to satisfy a 
c?-dimensional nonlinear Schrodinger equation (NLS) if 

idtu = —An + k\u\ p ^ 1 u, 

where A is the Laplacian operator in p is the nonlinearity parameter, 
and k is a parameter which is either +1 or —1. The case of interest in this 
article is p = 3 and k = — 1, called the 'focusing cubic NLS'. The focusing 
cubic NLS is an equation of interest in nonlinear optics, condensed matter 
physics and a number of other areas [I5|infll22|H0|l2|II7|l551l3S]. 

The questions of local and global well-posedness of nonlinear Schrodinger 
equations are still not completely understood. Local existence results un- 
der restrictive conditions on the initial data have been established in low 
dimensions |18| [21] , while ill-posedness results are known in higher dimen- 
sions |13j . For a survey of the literature and further references, see the recent 
monograph of Raphael |32j . 

One approach to this problem is via the method of invariant measures 
for the NLS flow, initiated by Lebowitz, Rose and Speer [21] and developed 
by McKean and Vaninsky [261 EH EH] and Bourgain j3[ 0J EJ [6] . Invariant 
measures coupled with Bourgain's development of the so-called X S)h spaces 
('Bourgain spaces') for constructing global solutions has led to important 
advances in this field. A striking recent development is the work of Tzvetkov 
[37] who used invariant measures and Bourgain's method to construct global 
solutions of certain nonlinear Schrodinger equations for random initial data. 
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The technique was further developed by Burq and Tzvetkov [TUl E] and Oh 
et. al. P31ET]. 

Most of the above works use invariant measures as a tool for proving lo- 
cal or global well-posedness for various classes of initial data. However, the 
nature of the invariant measures themselves have not been so well studied. 
Studying the nature of the invariant measures may yield important infor- 
mation about the long-term behavior of these systems. The only results 
we know in this direction are those of Brydges and Slade [7j in d = 2 and 
Rider |33t 151] in d = 1 . Some progress for invariant measures of the KdV 
equation has been made recently in [30]. In this article we investigate the 
case of the discrete NLS in three and higher dimensions. Two problems that 
immediately arise are: (a) the construction of the invariant measure for the 
cubic NLS due to Lebowitz, Rose and Speer [21] does not give a meaningful 
probability distribution when d > 3, and (b) local and global well-posedness 
are not well-understood in d > 3. To overcome these (very difficult) hurdles, 
we drastically simplify the situation by discretizing space and considering 
the so-called discrete NLS. Several well-posedness results under general con- 
ditions on the initial data are known for the discrete system [36, 38J, and 
the existence of the natural invariant Gibbs measure is straightforward. 

In return for this simplification of the problem, we give a large amount of 
refined information about the nature of the invariant measure. In particular: 

• We 'solve' the system 'exactly' in the sense of statistical mechanics 
by computing the limit of the log partition function. 

• Analysis of the partition function yields a first-order phase transi- 
tion; we identify the exact point of transition. 

• We prove the existence of so-called localized modes (also called dis- 
crete breathers |1 T|, \25\ 138] ) in functions drawn from the invariant 
measures and compute the size of these modes. 

• Additionally, we show that the localized mode persists at one site 
for an exponentially long time. 

The results are stated in Section [2] and the proofs are presented in later 
sections. The proofs involve elementary probabilistic arguments. 

Incidentally, we do not know how to take our results for the discrete NLS 
to some kind of a 'continuum limit' as the grid size goes to zero. Neither 
do we know how to derive conclusions about the long-term behavior of solu- 
tions of the discrete NLS from our results about the nature of the invariant 
measures. As of now, these are open problems that may well be unsolvable. 

2. Set-up and statements of main results 

It is of interest to study the nonlinear Schrodinger equation on graphs 
[U [H], and this is the setting of our results. 

Let G = (V, E) be a finite undirected graph without self-loops. Let D be 
the maximum degree of G and let n = \ V\ be the size of the graph. Let h 
be a positive real number, denoting the 'distance' between two neighboring 
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vertices in G. When G is a part of a lattice, h is called the lattice spacing 
(e.g., in [38])- 

For example, when G is a discrete approximation of the <i-dimensional 
unit torus [0, l] d represented as {0, 1/L, 2/L, . . . , (L — 1)/L} d , then n = L d 
and h = 1/L. In particular, n — > oo as L — > oo; moreover, if d > 3, nh 2 also 
tends to oo as L — > oo. This last condition will be crucial for us. There is 
nothing special about the torus. The condition nh 2 — > oo should hold for 
any nice enough compact manifold of dimension > 3. The condition is also 
satisfied if instead of the torus we take a cube whose width is increasing to 
infinity as the grid size 1/L tends to zero (since h = 1/L and n > L rf in 
this scenario), which can be viewed as a discrete approximation to the whole 
of R d . 

The condition nh 2 — > oo is satisfied in 2D and ID only if the domain that 
is approximated also tends to an infinite size. For example in ID, if the grid 
size is h = 1/L for some large L (tending to infinity) then the condition 
nh 2 — > oo holds on an interval [— K, K] only if K grows to infinity with L 
so fast that K(L)/L — > oo. 

The discrete nearest-neighbor Laplacian on G is defined as 



h 2 

y^x 

where y ~ x denotes the sum over all neighbors of x and / = (f x )xev is an Y 
map from V into C. Note that the scaling by h 2 is meant to ensure, at least 
in the case of the d-dimensional torus, that the discrete Laplacian converges 
to the true Laplacian as the grid size goes to zero. 

The discrete cubic NLS on G with the discrete nearest-neighbor Laplacian 
A is a family of coupled ODEs with f x = f x (t) : 

i-jifx = -A/x + n\f x \ 2 fx, x £ V, 

where k may be +1 or —1. These are known as the focusing and defocusing 
equations, respectively. The object of interest in this paper, as in Lebowitz- 
Rose-Speer [23] and Bourgain [3J, is the focusing NLS: 

(1) ij t fx = "A/, - \f x \ 2 f x . 

The discrete Hamiltonian associated with the focusing NLS is: 

(2) H{f):= 2 - £ L J '' 2 ' Y. I-''- 



n 

{x,y)eE 



n 

x&V 



Up to scaling by a constant, this is simply the discrete analog of the contin- 
uous Hamiltonian considered in [24, |3| |26]. 

Let the power be N(f) := Ylxev l/z| 2 ; the mass is n~ 1 N(f), but we will 
use these terms interchangeably. Then we have conservation of mass under 



1 



SOURAV CHATTERJEE AND KAY KIRKPATRICK 



the dynamics, using standard finite-difference techniques 
as weli as conservation of energy: 

Hence by the Liouville theorem, the measure dfi := e~^ H ^ Ylx^v ^ x * s 
invariant under the dynamics of the discrete NLS ([I]) for any real j3. How- 
ever, this measure has infinite mass if /3 > 0. The problem is easily solved 
by a mass cut-off as in |24} [3] (allowed due to conservation of mass) and 
normalization. The resulting probability measure 

(3) dfi := Z~ l e-^h [N{mBn} [] df x 

continues to be invariant under the NLS dynamics. Here B is an arbitrary 
positive cutoff, and Z is the normalizing constant (partition function). Of 
course, both fi and Z depend on the pair (/3, B). Let ip be a random element 
of with law fi. That is, ij) is a random function on V such that for each 

F(iPcA) = Z- 1 [ e-P H Vh mf) < Bn} df, 

J A 

where df = Y[ x df x denotes the Lebesgue differential element on C v . Our 
objective is to understand the behavior of the random map ip. The first step 
is to understand the partition function Z. The first theorem below shows 
that if we have a sequence of graphs with n and nh 2 both tending to infinity, 
the limit of n" 1 logZ can be exactly computed for any positive /3 and B. 

The result can be roughly stated as follows. Let m : [2, oo) — > M. be the 
function 

10 F 2 , (I 1 L 2 

It may be easily verified that m is strictly increasing in [2,oo), m(2) < 
and m(3) > 0. Thus, m has a unique real zero that we call C . Numerically, 
B r ~ 2.455407. Let 



(5) F(J3,B) 



| log(Svre) if pB 2 < 9 C , 

1 log(57re) + m(f3B 2 ) if pB 2 > 6 C . 

(Figure [l] shows a graph of F versus /3 when B = 1.) Theorem 2.1 below 
asserts that if /3 > and we have a sequence of graphs such that n — > oo and 
nh 2 — > oo (as in the d-dimensional torus for d > 3), all other parameters 
remaining fixed, then 

lim =F(P,B). 
The theorem also gives an explicit rate of convergence. 
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P 

e' c =2.455... " 

Figure 1. The free energy is constant for small inverse tem- 
perature and starts increasing at the critical threshold. Here 
mass is normalized, B = 1. 



Theorem 2.1. Suppose j3 > 0. Take any e G (0,1/5). There exists a 
positive constant C depending only on e, f3, B, h and D such that if n > C, 
then 

_ B ) > -Cn- 1 ' 5 - C{nh 2 )- 1 

n 

and 

logZ _ jCn- 1 /^ + Cn~^ if flB 2 < 9 C , 

The behavior of the random map ■0 is the subject of the next theorem. It 
turns out that the behavior is quite different in the two regimes (3B 2 < 9 C and 
[3B 2 > 9 C . To roughly describe this phase transition, let M\{ip) and M2(V>) 
denote the largest and second largest components of the vector (\i/) x \ 2 )xev- 
It turns out that when f3B 2 > 9 C , there is high probability that Mi(ip) ps an 
and M2(tp) = o(n), where 



(6) a = ^^ : =f + fV 1 -^- 

In other words, when f3B 2 > 9 C , there is a single x where ip x takes an 
abnormally large value, and is relatively small at all other locations. More- 
over, N(ip) « Bn with high probability. A consequence is that the largest 
component carries more than half of the total mass: 

\^x\ 2 a 1 
max w — > -. 

x EJ^r B 2 

On the other hand, when /3B 2 < 9 C , then Mi(ip) = o(n), but still N(ip) « 
Bn. Consequently 
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Figure 2. The fraction of mass at the heaviest site jumps 
from roughly zero for small inverse temperature, to roughly 
.71 at the critical threshold. (Here B = 1.) 



(Figure [2] shows the graph of the fraction of mass a at the heaviest site versus 
j3 when B = 1.) When j3B 2 > 8 C , the energy density H(ip)/n is strictly 
negative and approximately equals —a 2 , whereas in the regime (3B 2 < 9 C , 
the energy density is close to zero. The formula for a shows that a does 
not tend to zero as f3B 2 approaches 6 C from above (in fact, it stays bigger 
than B/2), demonstrating a first-order phase transition. These results are 
detailed in the following theorem. 

Theorem 2.2. Suppose h = n~ p for some p G (0, 1/2). Let a = a(/3,B) be 
defined as in ^ and let Mi(tp) and M2(ip) be the largest and second largest 
components of (\ipx\ 2 )xeV • First, suppose (3B 2 > 9 C . Take any q such that 
max{2p, 4/5} < q < 1. Then there is a constant C depending only on (3, B, 
D, p and q such that if n > C , then with probability > 1 — e~ nq l c , 



(7) 



n 

Mi(V) 



n 



< C n-^l\ 



B 



n 



anc 



< Cn-^/ 2 , 



n 



Next, suppose /3B 2 < 6 C . Take any q satisfying max{(l + 2p)/2, 17/18} < 
q < 1. Then there is a constant C depending only on /3, B, D, p and q such 
that whenever n > C, with probability > 1 — e.- nq / c 



(8) 



■n 



< 2n 



-2(1-9) 



B 



n 



< n 



-(l-q) 



and 



n 



Finally, iff3B 2 = 9 C andq is any number satisfying max{(l+2p)/ '2, 17/18} < 
q < 1, then there is a constant C depending only on (3, B, D, p and q such 
that whenever n > C, with probability > 1 — e~ n<1 l c either ([7]) or Q holds. 
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An obvious shortcoming of Theorem 2.2 is that it does not give a precise 
description the critical case (3B 2 = 9 C . It is important to know whether (|7) or 
([8]) is more likely, and how much. Another deficiency of both Theorem |2.1 



and Theorem |2.2| is that they say nothing about the thermodynamics of 
the 2D cubic NLS. A substantial amount of information is known about 
the 2D continuous system (see e.g., [32] and references therein) but precise 



calculations along the lines of Theorems 2.1 and 2.2 would be desirable. 

Additionally, it would be nice to be able to extend the theory to other 
nonlinearities than cubic. 

Theorem 2.2 says that when [3B 2 > 6 C , there is a single site x £ V which 
bears a sizable fraction of the total mass of the random wavefunction tp. 



This fraction is nearly deterministic, given by the ratio a/B. Theorem 2.2 
also implies that this exceptional site bears nearly all of the energy of the 
system. This is because the total energy H(ip) is approximately —a 2 n, while 
the energy at x is, summing over just the neighbors y of x: 



n 



-n^M^) 2 + 0(h- 2 ) 



-a n + o(n) 



the equality by Theorem 2.2 Such a site is sometimes called a localized 
mode. 

It follows easily as a corollary of this theorem that typical discrete wave- 
functions above the critical threshold have divergent discrete H 1 norm: 

2 



2_ 
fl 1 



fx fy 



h 



n n ^-^ 

xGV (x,y)£E 

However, it is not so clear that the discrete H 1 norm diverges even if (3B 2 < 
9 C . The following theorem shows that the the divergence happens on the 
discrete torus in dimensions > 3 for all values of /3 and B. 



Theorem 2.3. Suppose the context of Theorem 2.2 holds, with h = n~ p for 
some p G (0, 1/2) and n > C. Let ifi = (?px)x£V be a discrete wavefunction 
picked randomly from the invariant probability measure ft defined in (|3j). // 
j3B 2 > 6 C , then there is a positive constant c depending only on (3, B, D 
and p such that 
hand, if f3B 2 < 



< cn p ) < e n whenever n > 1/c. On the other 
l c , then the same result holds with a small modification: 
< e~' 



\gi < c\/6ri p ) < e on , where 5 is the average vertex degree. 

Since the measure /i of (|3j) is invariant for the discrete NLS ([T]), one 
may expect from the above discussion that if the initial data comes from fl, 
localized modes will continue to exist as time progresses. The question is 
whether the same site continues to be a mode for a long time (in which case 
we have a 'standing' or 'stationary' wave with a localized mode, sometimes 
called a discrete breather), or not. The following theorem shows that indeed, 
the same site continues to be the localized mode for an exponentially long 
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period of time. This is an example of a dynamical result deduced from a 
theorem about the statistical equilibrium. Of course, we need to use the 
NLS equation ([!]) for some basic dynamical information at one point in the 
proof. (As a side note, let us mention that global-in-time solutions of ([I]) 
are known to exist 



Theorem 2.4. Suppose h = n~ p for some p € (0, 1/2) and j3B 2 > C . Let 
a be defined as in Let ip(t) = (ip x {t)) X £V be a discrete wavefunction 
evolving according to ([!]), where the initial data ip(0) is picked randomly 
from the invariant probability measure fi defined in Choose any q such 
that max{2p, 4/5} < q < 1. Then there is a constant C depending only on 
f3, B, D, p and q such that if n > C , then with probability > 1 — e _n< ' C 
the inequalities ([7]) hold for ip(t) for all < t < e nq l c , and moreover there 
is a single x € V such that the maximum of \ip y {t)\ is attained at y = x for 
all <t < e n<1 l c . In particular, ifj(t) is approximately a standing wave with 
localized mode at x for an exponentially long time. 

The above theorem proves, in particular, the existence and typicality of 
solutions of ([I]) that have unique stable localized modes for exponentially 
long times if the initial energy or mass are above a threshold. One key 
difference between this theorem and earlier results about existence of discrete 
breathers (e.g., [38]) is that the earlier results could prove the existence of 
localized modes only if the mass was very large, i.e. tending to infinity, while 



Theorem 2.4 proves it under finite mass and energy. Another difference 
is that it shows the typicality, rather than mere existence, of a breather 
solution. 

Our final theorem investigates the probability distribution of the indi- 
vidual coordinates of a random map V picked from the measure p,p t B- It 
turns out that it's possible to give a rather precise description of the dis- 
tribution for small collections of coordinates. If fiB 2 < 9 C , then for any 
xi,...,X). G V, under a certain symmetry assumption on G, the joint distri- 
bution of B~ 1 / 2 (ip Xl , . . . ,^x k ) is approximately that of a standard complex 
Gaussian vector, provided k is sufficiently small compared to n. When 
/3B 2 > 6 C , the same result holds, but for the vector (B — a) 1 ^ 2 (^ Xl , . . . , ip Xk ) 
where a is defined in 

The symmetry assumption on G is as follows. Assume that there exists 
a group S of automorphisms of G such that: 

1. |E| = n. 

2. No element of S except the identity has any fixed point. 

When these conditions hold, we say that G is translatable by the group S. 
For example, the discrete torus is translatable by the group of translations. 
Note that a translatable graph is necessarily transitive. 

Theorem 2.5. Suppose the graph G is translatable by some group of au- 
tomorphisms, according to the above definition. Suppose h = n~ p for some 
p £ (0,1/2), and let ip be a random wavefunction picked according to the 
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measure p,. Take any k distinct points x±, . . . , xt G V. Let <p = ((pi, . . . , 
be a vector of i.i.d. standard complex Gaussian random variables. If (3B 2 < 
9 C , then there is a constant C > depending only on (3, B, D and p such 
that if n > C , then for all Borel sets U C C k , 

\F(B~ l / 2 (ij; Xl ,. . . ,il>x„) € U) - F(<f> € U)\ < kn' 1 ' '. 

If f3B 2 > 9 C , the result holds after B^ 1 ^ 2 is replaced with (B — a)~ l l 2 where 
a = a(/3,B) is defined in and the error bound is changed to k?n~ 1 l c . 

Obviously, it would be nice to have a similar result for f3B 2 = 6 C , but our 
current methods do not yield such a result. 

When (3B 2 > 9 C , the reader may wonder how can ijj x behave like a com- 



plex Gaussian variable with second moment B — a, when Theorem 2.2 and 
symmetry among the coordinates seem to imply E|^ a .| 2 ~ B. It is exactly 
the peaked nature of the field in the case (3B 2 > 6 C which allows for a 
convergence in law without that of the second moment. 

The rest of the paper is devoted to the proofs of the theorems of this 
section. Some preliminary lemmas are proved in Section [3| Theorem |2.1| is 



proved in Section [4j Theorem 2.2 in Section [5j Theorem |2.3| in Section [6j 



Theorem 2.4 in Section [7j and Theorem 2.5 in Section [8 



3. Preliminary lemmas 

Throughout the rest of this article, C will denote any positive function 
of (e, (3,B,h,D), whose explicit form is suppressed for the sake of brevity. 
The value of C may change from line to line, and usually C will be called a 
'constant' instead of a function. When h = n~ p , C will depend on p instead 
of h. 

Some further conventions: Sums without delimiters will stand for sums 
over all x £ V. For each f £ C v and each k > 2, let S k (f) := E\fx\ k - In 
particular, N(f) = 82(f)- Define 

Z':= [ e^fh {S2(mBn} df 

If S 2 (f) < Bn, then 

0< \f*-fy\ 2 < E WA 2 + 2\f y \ 2 )<ABDn. 

(x,y)eE (x,y)&E 

Thus, we have the important bounds 

(9) e~ Ch ~ 2 Z' < Z< Z' 

which precipitate the irrelevance of the kinetic term in the Hamiltonian 
(although this is not obvious a priori). For each < a < 6, define 

T a ,b ■= if£C v : a 2 n 2 (l - n^) < S 4 (f) < aV(l + n~V 5 ), 



bn(l - n" 1/5 ) < S 2 (f) < bn(l + n" 1 / 5 )}. 
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Define the function 

(10) L(a, b) := (3a 2 + log(6 - a) + log vr + 1. 



It will be shown in Lemma 3.6 that in fact, F(/3,B) = sup < a<b < B L(a,b). 
Lemma 3.1. For any < a < b < B/(l + n~ 1//5 ) ; 

( r, ^ Cn4/5 C \ 

Z > a exp nLla, b) — ; —k . 

\ b — a n z J 

Proof. Fix < a < b < B/(l + n^ 1 / 5 ). Let = (^) xg V be a collection of 
i.i.d. complex Gaussian random variables with probability density function 



1 



. e -|*| 2 /(6-o). 



7r(6 — a) 

Note that (f> can also be viewed as a random mapping from V into C. Now, 

r l s 2 (/) 

Consequently, 

(11) Vol(T 0)6 ) > exp (^(1 - n" 1 / 5 ) + nlogvr(6 - o))p(^ G r a>6 ). 

Fix an element o G V. Define the sets 

£ : := {/ : on(l - n" 1 / 4 ) < |/ | 2 < an(l + n" 1 / 4 )}, 
^^{/rmaxl/.^aV/ 4 }, 

£ 3 := {/ : (b - a)n(l - n^ 4 ) < ]T | | 2 < (6 - a)n(l + n" 1 / 4 )}. 

Suppose / G -Ei Pi £2 H £3. Since / G £1 and / G £3, 

6n(l - n" 1 / 4 ) < S 2 (/) < 6n(l + n~ 1/4 )- 
Again since / G £1 and / G £2, if n > C, 

S±(f) < a 2 n 2 (l + n" 1 / 4 ) 2 + a 2 n 7 / 4 < a 2 n 2 (l + n" 1 / 5 ), 
and similarly 

5 4 (/) > a 2 n 2 (l -n- 1 / 5 ). 

Thus if n > C, then £1 n £ 2 n £3 C r a>6 . 

For any x £ V, the real and imaginary parts of 4> x are i.i.d. Gaussian with 
mean zero and variance (b — a)/2. Thus, \4> x \ 2 is an exponential random 
variable with mean b — a. The inequality e~ u — e~ v > (v — u)e~ v that holds 
for v > u gives 

/ an(l + n ' ) 

G £1 > -t exp ^ '- 

b—a \ b—a 



PROBABILISTIC METHODS FOR NLS 11 

Further, note that by a simple union bound 

F(<P # E 2 ) < {n - l) e -°« 3/8 /(6-«) 
and by Chebychev's inequality, 

Thus, if n > C, then F(<p G E 2 n E 3 ) > 1/2. Lastly, observe that the event 
{4> G Ei} is independent of {0 G i?2 H -E3}. Combining these observations, 
we see that when n > C, 



g r 0j6 ) > P(0 g Bi)P(0 e£ 2 n # 3 ) 

(12) / an Cn 3 / 4 x 

> aexp - . 

V — a — a / 



This inequality, together with gives that if n > C, 

Vol(r a>6 ) > a(evr(6 - a)) n e" Cn4/5/(6 - a) . 
Since 6 < + n" 1 ^), r a>b C {/ : 5 2 (/) < Bn}. Thus, if n > C, 

Z'> [ e^~ ls ^df 

> e^ 2 "( 1 -" 1/6 )voi(r 0i6 ) 

> a(e 1+/3a2 vr(6 - a)) n e- Cni/5 = ae nHa,byCn^/(b-a) 

This lower bound and Q complete the proof of the lemma. □ 
Lemma 3.2. For any ra~ 2e / 5 < a < b < B , 

C n 4/5+e 



Z' ab := [ e pn 1S4 ^df < exp( nL(a, b) + 



b — a 



Proof. Fix n 2<E / 5 < a < b < B. Let 4> be the random vector defined in the 
proof of Lemma |3.1| Let 



T 'a,b ■= U {/ G r ^ : on (! " Su" 1 /^) < |/ | 2 < on(l + n" 1 / 5 ^), 
max I / a 1 2 < 4an 4 / 5+e ). 



Recall that for any o, |^ | 2 is an exponential random variable with mean 
b — a. Thus, 



Gr^<£P(|^ | 2 >an(l-3n-V5+<)) 



(13) ° eV 

/ an Cn A / 5+e 

= nexp ( - h —r 

V b — a b — a 



A 2 := {/ : 317 C V, \U\ < n^ +2 ^l\ E \f x \ 2 > an(l + n 
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Next, define 

M. :={/:|{^:|/ I |>n( 1 - £ )/ 5 }|» i ( 4+2£ »/ 5 }, 

-l/5+e- 

We claim that if n > a _5 ' 2e (which is true by assumption), then 

(14) r^cr^uijuia, 

To see this, take any f £ A^n A^d T a ^. Let 

[/:= {x:\f x \>n^ 5 }. 
Since / € \U\ < n (i+2e ^ 5 . Therefore, since f € A%, 

(15) ^IM^a^l + n- 1 / 5 ^). 

set/ 

Again since / 6 ^ we a l so have 

(16) max|/ a .| 2 < an(l + n" 1/5+e ). 

Note that since n > a~ 5 / 2e , 

E l/*| 4 = E IM* ^ n 1 ^ 1 -)/ 5 < a 2 n 9 / 5 . 

Since / G r ajb , 

ElM 4 >aV(l-n-V5). 

Thus, 

> a 2 n 2 (l - n-V5 _ „-i/5) = a 2 n 2 (1 _ 2n -i/ 5) _ 
On the other hand by ( |15[ ), 

El/xl^^axl/^nElM 2 

< (maxl/^^a^l + n- 1 / 5 ^). 

Combining the last two displays implies that 

, , l2 a 2 n 2 (l - 2n- 1 l r °) 

max \j x \ > — , 

i an(l + n i / 5+e ) 

> on(l - 3ra- 1/5+e ). 



Together with (16), this shows that 

(17) on(l - Sn" 1 / 5 ^) < max \f x \ 2 < an(l + n" 1/5+e ). 
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Next, let o be a vertex at which / attains its maximum modulus. Let x 
be any other vertex. If x ^ U, then since n > a~ 5 / 2e , |/ x | 2 < n 2 *- 1- ^/ 5 < 
an 2 / 5 < an 4 / 5+e . If x G £7, then by (pi), 



|/o| 2 + |/,| 2 <^|/ y | 2 <an(l + n- 1 / 5+e ), 



and therefore by (17), 



\f x \ 2 < an(l + n" 1 / 5 ^) - an{\ - Sn" 1 / 5 ^) = 4an 4 / 5H 



Hence /GT^ from the above display and (17), and the claim (14) follows. 
Consequently, if n > a~ 5 ^ 2e , 



(18) 



6 r a , 6 ) < P(0 g r' 6 ) + P(0 g A x ) + ¥{</> g A 2 ). 



If G Ai, then there is a set C/ C V such that |J7| = [n (4+2e) / 5 ] and 
|<^a;| > n^ 1- ^/ 5 for all x £ U. Therefore, a union bound over all possible U 
gives 



(19) 



e j,U| n )f e -» 2(1 - i)/5 /(M 
lj " 1 [n( 4 + 2 ^)/5] 1 1 



n (4+2e)/5 



< exp logn 



n 



6/5 



b — a 



Take any [/ C V and let j := |C/|. Then (6 — a) 1 X^ef/ 1^1 2 * s the sum of 
j i.i.d. exponential random variables with mean 1. Thus, for any t > 2, 



X^l 2 > (b~a)t 



(7-1)1 



e 31 cte 



(x + t) 



3-1 



-e x ote 



(J " 1)! 

/•°o 2 ^ 2 (^- 1 +^- 1 ) 

< e - t (2 J '- 2 + ^'- 1 ) <tV"*. 



On the other hand, for < t < 2, the bound Ce * works. If <j) £ ^2 then 
there exists U QV with = [n( 4+2e )/ 5 J and J2 x &u \f*\ 2 ^ ^(l+n" 1 /^^. 
Thus, with t = (b - ay 1 an{l + n" 1 / 5 ^) and j = [n (4+2e)/5 J , the above 
inequality, the assumption that a > n~ 2e / 5 , and a union bound over all 
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possible U shows that if n > C, 



n 



eA 2 )< I „ f4+2 , /RI h— : ex P 



(4+2e)/5 -, /t . , 



j an(] 



[n^)^\)\b-aj F V b-a 

< ewfcnW log - an(1 t n " V5+£) l 
\ b — a b — a J 



n an n ( 4 + 3e )/ 5 



< exp 

\ b — a b — a b — a 

, . ( n( 4+3e )/ 5 an 

\ 2(b — a) b — a / 

Together with ([13]), ([18]) and (fT9|), this shows that if n > C (and n > a~ 5//2e ), 
then 

fCn^ 5+e an 
e r a)6 J < exp - 



V 6 — a 6 — a 



Now, similar to the beginning of Lemma 3.1 



1 

7r n (6- a) 



s 2 (/) 
e b ~ a (if 



^ 1 -5t^(1+™ _1/5 )a^ i/r ^ 

^ — w e b ~ a ; Voi(r o6 ) 



Therefore, if n > a 5//2<E , 

Vol(r 0j6 ) <exp(- 5r ^(l + n- 1 / 5 ) + nlog^(6-a))p(0er aife ) 
(21) \b-a J 

<(evr(6-a))^ Cn4/5+E /(^). 
Finally, observe that 

Z' ab = [ e^- ls ^df < e^ 2 "( 1+ "" 1/5 )Vol(r a b ). 

This completes the proof of the lemma. □ 
Lemma 3.3. For each n > 1 and r > 0, let 

U n , r := {/ € C" : S 2 (f) < rn}. 
Then Vol(l7 ft , r ) < (rve) n = e nL( -°> r \ 
Proof. Note that 

1 = / r ^ d f 

f e S 2 (f)/r i 

> / 7 , ra df > -Vol(C/ n , r ). 

(The last inequality holds because 82(f) < ^n on f7„ )? ..) □ 
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Lemma 3.4. For each f G C v , let a{f) := y/S^J)/n and b{f) := S 2 (f)/n. 
Note that a <b. For each c > let 

W c := {/ e C v : &(/) < B, b(f) - a(f) < c}. 

Then 



f e Pn- l S 4 (f) d f < e (C+log C )n^ 

JWc 



Proof. Take any / 6 W c . Let a = a(f) and 6 = b(f). Then 

max | f |2 > T, x \fx\ A _ o^n 

Thus, if o is a vertex at which |/| is maximized, then 
]T|/,| 2 < ^-^n< 2(6- a)n. 

Let d := cn/(n — 1). The above inequality and Lemma 3.3 show that 
Vol(W c ) < £ Vol({/ : \f Q \ 2 < Bn,Y,\f*\ 2 < 2cnX) 

= B7rn 2 Vol(C/ n _ li2c 
< Svrn 2 (2c / vre) n - 1 . 

Since 

/ e? n ~ ls ^df < eP B2n Vol(W c ), 
Jw c 

this completes the proof. □ 
Lemma 3.5. For each a > 0, /ei 

V a ■= {/ : &(/) < B, a(f) < a}, 

Then 

I e/ 3n- 1 5 4 (/) (j y < e /3a 2 n+L(0,B)n^ 

JVa 

Proof. Simply note that 

e ^n- 1 5 4 (/) d y < e ^ 2 "Vol([/„ iB ), 



and apply Lemma 3.3 □ 

Lemma 3.6. For each 9 > 0, let fg : [0, 1) — > M. be the function 

f e {x) := 6x 2 +\og{l-x). 

Then there is a 9 C > such that if 6 < 9 C , fg has a unique maximum at 
x = 0, whereas if 9 > 9 C , then fg has a unique maximum at the point 
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When = 6 C , fg attains its maximum at two points, one at zero, and the 
other x*(9 c ). The number 9 C is the unique real solution of 



9 19 2 , 1 1 

i-j + aV 1 -« +1 '*(2-aV 1 -» / J- a 

Numerically, 9 C ~ 2.455407. Lastly, ifF((3, B) is the constant defined in 
then 

F(j3,B)= sup L(a,b) = sup L(a,B), 

0<a<b<B 0<a<B 

and 

(22) L(a, B) = log(S7re) + f pB> (a/B). 

Proof. Let Mq := sup 0<;c<1 fe(x). It is easy to see that Mq = 0, and Mq > 
for sufficiently large 9. Moreover Mq is an non-decreasing function of 9. Let 

9 C := inf{# : Mq > 0} = sup{# : Mq = 0}. 

If 9 < 9 C , then we claim that fg has a unique maximum at zero. It is 
clear that zero is a point of maximum. To show that it is unique, suppose 
not. Then there is some x > where fe(x) = 0. Then for & 9' £ (9,9 C ), 
fg'(x) > 0, giving a contradiction to the definition of 9 C . 

When 9 > 9 C , it is clear that the maximum must be attained at a non-zero 
point. To show that it is unique, observe that 

f>( x ) = 29x-- i —, f' e '(x) = 29 ' 



1-x J ° y ' (1-x 



,2 ' 



and therefore there may be at most two points in (0, 1) where fg vanishes, 
and exactly one of them can be a maximum. Solving the quadratic equation 
shows that the maximum is attained at x*(9) (and also shows that 9 C > 2). 
From the formula, it is clear that when 9 — > 9 C , the limit of x*(9) exists and 
is positive. By uniform convergence on compact intervals, it follows that 
fe c (x*(9 c )) = 0. The equation fe c (x*(9 c )) = is equivalent to the defining 
equation for 9 C . It is a unique solution of the equation because the left hand 
side can be easily shown to be a strictly increasing function of 9 in [2,oo). 
Finally, note that for any < a < b < B, L(a, b) < L(a, B), and the easy 



identity (22) implies the relation between F and L because m(9) defined 
in Q is nothing but M e if 9 > 9 C . □ 



4. Proof of Theorem 2.1 (Partition Function) 



Let be defined as in Lemma 3.6. Recall the equation (|22|) , which implies 



that L(a, B) is maximized at Bx*(a/ B). Let us call this point a* . (The point 



a* is called a(/3, B) in the statement of Theorem 2.2 ) Note that L(a, b) is 



maximized at (a* , B), since L(a, b) < L(a, B) for any < a < b < B. 



By Lemma 3.6, we know that a* = if f3B < 9 C and a* > when 
f3B 2 > 9 C . Moreover, when f3B 2 ^ 9 C , a* is the unique point of maximum. 
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4.1. Supercritical Partition Function. First, consider the case f3B 2 > 
9 C , so that a* > 0. By Lemma 3.1 it follows that if n > C, 

^ > L(a*, + n-VB)) - Cn- 1 ^ - 

But L has bounded derivatives in a neighborhood of (a*,B). Thus, if n > C, 

(23) ^ > L{a\B) - Cn- 1 / 5 - ~. 

Since L(a*, B) > L(0, £?), there exists ao = ao(/S, 5) small enough such that 

,2 



where 



Pa z + L(0,B) < L(a*,B)-5, 



6:=-(L(a*,B)-L(0,B)). 



Therefore by Lemma 3.5 

(24) f e P™- 1 S4f) d f < e n(L(a*,B)-S)^ 

JVa 

Since logc — > — oo as c — > 0, there exists Co = Q)(/3, S) > such that 

C + logco < L(a*,B) - 5, 



where C is the constant in Lemma 3.4 Therefore by Lemma 3.4 



(25) 



Wc 



Now take any / such that b(f) < B. (Recall the definitions of a(f) and b(f) 
from the statement of Lemma |3.4| ) Let A be a finite collection of (a, 6) such 
that 

|J T a , 6 D {/ : b(f) < B, a(f) > a , &(/) " «(/) > c }- 

(a,b)e.4 

It is easy to see that if n > C, then „4 can be chosen such that |.4| < Cn 2 / 5 , 
and for each (a, b) £ A, a > ao/2 and b — a > cq/2. Let us choose such a 
collection .4,. Then by Lemma 3.2, for each (a, b) £ .4, 



7' <- „nL(a,b)+Cn 



4/5+6 



(26) 

Now, if / and / Wco> then a(_f) > ao and b(f) — a(f) > c$. Therefore 
by ([24]), @ and po]), it follows that if n > C, 

Z' < 



= /3n- 1 S 4 (/) 



V ao UW co 



df+ E 



(o,6)£v4 



< 2e "( i ( a *' B )- 5 ) + |^| max e n(2*«,&)+Cn-V«+«) 

(a,b)g.A 



< 



Cn 2/5 e n(L(a* ,B)+Cn^/ 5 +^) _ 



A combination of the above inequality, (23), ([9]) and Lemma 3.6 proves the 
conclusion of Theorem 2.1 when f3B 2 > # c . 
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4.2. Subcritical Partition Function. Next, consider the case fiB 2 < 9 C , 
when L(a, b) is maximized at (0, B). (The point of maximum is not unique 
when (3B 2 = 9 C , but that will no t ma tter in the proof.) Taking a = n -1 / 5 
and b = B/(l + n _1//5 ) in Lemma 3.1 it follows that 



^ > Lin' 1 / 5 , B/(l + n" 1 / 5 )) - Cn 1 ' 5 
n 



C __ Clogn 



n 



n 



Since L has bounded derivatives in a neighborhood of (0, B), this shows that 
if n > C, 



(27) 



^lyL^B)-^- 1 / 5 
n 



C_ 
h 2 ' 



n 



By Lemma 3.4 a constant c\ = ci(e, /3, B) > can be chosen so small that 
when n > C, 



(28) 



Lemma 3.2 shows that whenever n 2e / 5 < a < b < B and b — a > c\, 

(29) l ° gZ a,b < L ^ 6) + Cn -l/5+e < L(Q) ^ + Cn - 



-1/5+e 



Let ai := n 2<E//5 . Let „4. be a set of (a, 6) such that 

IJ T a , b D {/ : &(/) < 5, &(/) - a(f) > Cl , o(/) > 2oi}. 
(a,&)e,4 

From the definition of r 0) &, it should be clear that A can be chosen such 
that \A\ < n c , and for each (a, b) *E A, a > ai and 6 — a > c\/2. Therefore, 



(29) holds for every (a, b) E .4. Thus, from (28), (29) and Lemma 3.5 we get 



< e 4/3a 2 in +L(0,B)n + g iZ(0,B)n + |^| e L(0,B)n+Cn 4 / 5+e 

Since |^4| < n*^ and a 2 n = n 1_4t//5 , this completes the proof in the case 
I3B 2 <9 C . 



5. Proof of Theorem 2.2 (Gibbs Measure) 



Recall that h = n- p in this theorem, where p E (0, 1/2). In this proof, 
whenever we say "for all (a, b) satisfying..." it will mean "for all (a, b) such 
that < a < b < B, satisfying...". 
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5.1. Supercritical Gibbs Measure. First, consider the case f3B 2 > 9 C . 
Choose q such that max{2p, 4/5} < q < 1. Choose e satisfying 4/5 + e = g. 
Note that e 6 (0, 1/5), as required. Let r := (1 + q)/2. Let a*, ao, Co, S, 
a(f) and &(/) be as in the proof of Theorem |2.1| Let 

A 1 := {/ : b(f) < B — 2n~^} D V a c n T^ c c , 
^2 := {/ : |a(/) - a*\ > 2n~^' 2 } n F a c n ^ c c . 

Let .Ai be a collection of (a, 6) such that < a < b < B and 

Clearly, A\ can be chosen such that \A\\ < Cn 2//5 , and for all (a, 6) E Ai, 
a > ao/2, b-a> co/2 and either 6 < B-n~( l ~ r ^ or |a — a*| > n~( 1 ~ r )/ 2 (or 
both). With these properties of Ai, two conclusions can be drawn. First, 
by Lemma 3.2, one can conclude that for each (a, b) G Ai, 



(30) Z' a b < e nL (^ b )+Cn 



4/5+e 



Next, observe that by equation (22) and Lemma 3.6, at the point (a, b) = 
(a*,B), dL/da = and d 2 L/da 2 < 0. Also observe that at all points (a, B) 
such that B — a > cq/2, dL/db is uniformly bounded away from zero and 
positive. Moreover the derivatives of L are continuous, L is increasing in b, 
and L(a,b) is uniquely maximized at (a*,B). 

Combining all of this, it follows that if n > C, then for all (a, b) such that 
\a-a*\> n-( 1 " r )/ 2 and b - a > c /2, 

n -(i->") 

L(a,6) < L(a,B) < L(a* , B) -— , 

and for all (a, b) such that b < B — n~( 1_r ) and b — a > cq/2, 

„-(l-r) „-(l-r) 
L(a, 6) < L(a, B) — < L(a* , B) — . 

Thus, if n > C, then for all (a, 6) G Al, 

L(o,6) < L(a*,B) - 



C 

From the above display, inequality ( |30[ ), and the fact that 4/5 + e < r, it 
follows that if n > C, then for all (a, 6) G Ai, 

(31) ^ 6 < e r^(«*,-B)-« r /C. 
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Let A := K U W Co UiiU A 2 . Then by (|31j), (|24j), (|25j), (|23j), (|9|) and the 
observation that r > max{4/5 + e, 2p}, we get that if n > C, 

P(V eA)<Z- 1 ( [ e Pn ~ ls ^df + \Ai\ max ^ b ) 

\Jv ao UW co (a,6)6^i ' / 

< e -nL(a*,B)+Cn 4 / 5 +Cn 2 P^ e ™( i ( £t *. B )-' 5 ) + <7 n 2/5gnL(a*,B)-n r /C^ 

Thus, with A' := {/ : b(f) <B, f & A}, we have that for n> C, 

(32) F(%b eA')>l- e~ nr/c . 
Note that if / G A', then 

(33) K X iV(/) - B| = - B| < 2n-( 1 - r ) 
and 

|n _1 fr(/) + a* 2 \ < Cn"* 1- ^ + \n~ 2 S 4 (f) - a* 2 \ 

<Cn-^ + \a(f)-a*\\a(f) + a*\ 

{:) '' ! < Cn-^- 2 '^ +2B\a(f) - a*\ 

< Cn^l 2 . 

Let ^.2 be a collection of (a, b) such that 

(35) |J r ai6 D^. 

From the definition of A' , it should be clear that if n > C, then .4.2 can 
be chosen such that < Cn 2 / 5 , and for each (a, 6) G ^.2, b — a > co/2, 
a > ao/2, 

(36) <4n" (1 ^ r,) , 
and 

(37) |o-o*| <4n-( 1 " r J/ 2 

Take any (a, 5) G -4.2- Let T' a b and be defined as in the proof of Lemma 
Then from Q, @ and @, it follows that if n > C, 

n 4/5+e 



3.2 



g r 0)6 \rl 6 ) < exp 



an 
b — a 



C 



(The exponent (4 + 3e)/5 in (20) improves to 4/5 + e since a > ao/2 here. 
This is easy to verify in the derivation of (20).) Therefore, as in (21), and 
using the fact that b — a > co/2, 



voi(r 0j6 \r' b ) < exp 



bn 



(1 + n" 1/5 ) +nlogvr( 



b — a 

< e n(L{a,b)-/3a 2 )-n 4 / 5 +' /C 



s r 0)6 \r 6 ) 
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Thus, by (23) and ^ and the observation that 4/5 + e = q > 2p, 
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P(^€T atb \T' atb )<Z- 1 J 



= /3n- 1 S 4 (/) 



df 



(38) 
Let 



<Z- l eP™ a ^- 1/6 )Vol(T a j b \T' atb ) 



Q== U r « 

(a,b)S^2 



6- 



By (35), 



^'^u( |J (r 0j6 \r' ai6 )). 

Therefore by the fact that \Ai\ < Cn 2 / 5 and d3§b, 



G A') < 



(39) 



eQ) + \A 2 \ max P(^r 4 , t \r; k ) 



< F(i/) G Q) + Cn 2/ V n 



.4/5+. 



By (39) and (32) we see that for n > C, 

G A') + P(V> G Q) — 1 

-n«/C 



-n«/C 



P(V> g A! n Q) > 

> 2P(V> G A') - e"^ /0 - 1 

> 1 _ e -nVC _ e -n'/C > J _ e 

Recall that M\{f) and M 2 (f) denote the largest and second largest compo- 
nents of the vector (\f x \ 2 )xeV- Note that if / G Q, then / G r' a b for some 
(a, 6) £ A 2 - The definition of r' a b implies that 

\Mi(f) - an\ < Cn^ 5+e , M 2 (f) < Cn 4/5+e . 
Combining this with ( |37[ ) and the fact that 4/5 + e = q gives 

|Afi(/)-o*n| < Cn (3+9)/4 . 



The last two displays and the inequalities ( |33| ) and (34) complete the proof 
of the theorem in the case f3B 2 > 6 C . 

5.2. Subcritical Gibbs Measure. Next, consider the case f3B 2 < 9 C . Let 
q satisfy 

max{(l + 2p)/2, 17/18} < q < 1. 

Choose e such that 

1 - — = q. 



Since g > 17/18, it follows that e < 1/7 and therefore 1 - 2e/5 > 4/5 + e, a 
fact that will be needed below. 

If M\(ij)) > n 1-26 / 5 , then S<i(ifi) > n 2_4<E//5 and therefore ip G T atb for s ome 
n" 2e / 5 < a < b < B. Let ci be as in (|28|) . Then combining Lemma |3.2 
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(27) and the facts that L is uniquely maximized at (0,B) and dL/da < 
at (0, B), we see that if n > C, then for any (a, b) satisfying 2n _2e / 5 < a < 
b< B and b — a > c%/2, 



m> e r a , 6 ) < 


z 


-i 7 t 




< 


e" 


nL(0,B)+C'n 4 ^+Cn 2 P- 


fnL(a,b)+Cn 4 / 5+E 


< 


e~ 


nL(0,B)+Cn 4 ^ 5 +Cn 2 P- 


fnL(a,B)+C*n 4 / 5+e 


< 


e" 


C -l n l-2 f /5 +Cn 4/5+ e4 





By our choice of e, it follows that for any such (a, b), if n > C then 

(40) P(V E r o , 6 ) < e-" 9 / c 

As usual, we can choose a set A3 of size < n c such that 

U ?a,b 5 {(a, 6) : 2n~ 2t/5 <a<b<B,b-a> c x }. 

Moreover, we can ensure that b — a > c\/2 and a > n~ 2e / 5 for all (a, b) E .A3. 
Therefore by (40), (27) and our choice of ci, we see that if n > C, then 

er ai6 )4-P(^w, 

< e -n"/C + z -X e \nL{Q,B) 

< e~ nq ' c . 



P(Mi(V>) > n 1_2e/5 ) < I A 3 1 max 

(a,6)e^ 3 



(41) 
Define the sets 



Ei 
E3 



if 

if 
{/ 



S 2 (f)>Bn-n«}, 
\H(f)\<2n 2 i- 1 }, 
Mx(f) < n*}. 



First, observe that if n > C, then £3 C E 2 : if M\{f) < n 1 2e / 5 = n q and 
n > C, then 

|ff(/)|< Cn^ + n^MiC/) 2 

< Cn 2p + n 2 « _1 

< 2n 2 *-\ 

since 2q — 1 > 2p. In particular, if n > C, 

E £ 2 ) = P(\H(ip)\ < 2n 2 «- 1 ) > 1 - e- n "l c . 
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Next, observe that by (27) and Lemma 3.3, the probability of ip belonging 
to El n E 2 is: 

P(S 2 0/0 <Bn- n q , \H(ij>)\ < 2n 2q ~ 1 ) 
= [ e-^fh {S2(mBn}d f 

J{f:S 2 (f)<Bn-nl, \H(f)\<2n 2 i-^} 

< Z-W^-'VoliU^Bn-rfi) 

< exp ( - n(L(0, B) - L(0, B - n -(1-ff) )) + Cn 4 / 5 + Cn 2p + Cn 29 " 1 

< e -ni /C+Cn^ 5 +Cn 2 f+Cn 2 i- 1 < g-n 9 /C 



Combining the last two displays with (41) gives that for n> C, 



e£infi 2 n E 3 ) = eB 2 n £ 3 ) - P(V> efi[n£ 2 n £ 3 ) 
= F(ij> e E 3 ) - e n£ 2 n £ 3 ) 

> 1 - e- n " /c - F(ip e E{ n £ 2 ) 

> 1 - e~ nq / c . 
This completes the proof in the case j3B 2 < 6 C . 

5.3. Critical Gibbs Measure. Finally, when j3B 2 = 9 C , L(a,b) is maxi- 
mized at exactly two points: (0, B) and (a*, B). As before, this implies that 
with high probability, ip cannot belong to any r a & where (a, b) is away from 
both of these optimal points. It can be deduced from this, exactly as before, 
that with high probability either ([7]) or ^ must hold. We omit the details. 
This completes the proof of Theorem \2.2\ 



6. Proof of Theorem 2.3 (Blow-up of H 1 norm 



The conclusion for the supercritical case f3B 2 > 9 C is almost immediate 
from 0. To see this, note that if x is the mode, then Y2 y ~ x — ^y\ 2 > 
an — o(n) with high probability. 

Consider the subcritical case j3B 2 < 9 C . Define E 2 := {/ : \H(f)\ < 
2n 2q ~ 1 }. Then by th e low er bound on Z from Theorem 2.1 and the subcrit- 

if n > C, then for any A C C , 



2.2 



ical part of Theorem 

eA)< F(if> E 2 ) + £ An E 2 ) 

e-^h {S2(mBn] df 



< e -n'/C + z -l f 

J A 

< e~ nq / c + (B7re)- n e 



AnE 2 

n„Cn 4 / 5 +Cn 2 P+C'n 2 i- 1 



AnE 2 



1 {S 2 {f)<Bn}df- 
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Observe that if is a random vector whose components are i.i.d. standard 
complex Gaussian, then for any U C C v , 



(B7re)- n [ l {S2 (f)<Bn}df < (^e)~" / e n ~ S * {f)/B l { s 2 (f)<Bn}df 
Ju Ju 

< (Bn)- n [ 
Ju 



e S2(f)/B df 
eU). 

Since the graph G has maximum degree D, it is easy to see (e.g., by a greedy 
algorithm) that there is a subset of edges E' C E such that \E'\ > \E\/D 
and any two edges in E' are vertex disjoint. Therefore, for a sufficiently 
large C (depending on the usual parameters) , using the fact that | <p x — 4> y \ 2 
is an exponential random variable with mean 2, 

^ e -Y. {x , y )eE\4>x-<t>y\ 2 ^ < E( e - \<t>x-<t>y?\ 

= \{ E(e-l*--^l a ) = 3-1^1 <e-\ E \' c . 

(x,y)eE> 

Thus, there is a constant C sufficiently large such that 

E \^-<f>y\ 2 <\E\/c)<e-W/ c . 

\x,y)eE ' 

Combining this with the two earlier observations and the observation that 
5 = 2\E\/n completes the proof in the case (3B 2 < 6 C . 

Finally, for the critical case, when f3B 2 = 9 C , either ([7]) or ^ must hold by 
Theorem 2.2 If ([7]) holds, the result is trivial. For functions satisfying ([8]), 
an argument similar to the above has to be made with small modifications. 
(For example, instead of a single set E%, we have to consider two sets Ei 
and E'2 such that E<i U E' 2 has high probability.) We omit the details. 



7. Proof of Theorem 2.4 (Standing Wave) 

As usual C denotes any constant that depends only on [3, B, D, p and 

so as 



q. Let Cq denote the constant C from the first part of Theorem 2.2 



not to confuse it with other C's in this proof. Let A denote the interval 
[Q,e n "/ 2Co ]. Define 

T := {t G A : V(i) violates @ for C = C }. 

Then by Theorem |2.2| and the invariance of the flow, if n > Co, 

E|T| = [ P(i € T)dt < e -n q /C Oe n«/2C = e -n"/2C ^ 
JteA 

where \T\ denote the Lebesgue measure of the set T. Therefore, 

P(|T| > e -«V4C 0) < _ni\ < e -nV4C . 
vl 1 — ' — p—ni/ACo — 
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Also by Theorem [Z2| P(0 G T) < e~ n9/c °. Thus, 

P(|T| < e -" 9 /4C an d T) > 1 - e "^/4C _ e -»«/Gi>. 

Suppose this event happens, that is, |T| < e - n9 /4Co anc i T. We 
claim that under this circumstance, the discrete wavefunction ip cannot have 
modes at two different locations for two distinct times t, s £ A and moreover 
ip(t) satisfies ([7]) for all t G A. This would complete the proof. 

By Q, we know that if t £ T, then for each y 6 1/, |Vv(*)| 2 must either 
belong to the interval [ore — Con^ 3 " 1 " 9 )/ 4 , are + Con^ 3+IJ ^ //4 ] (such a point will be 
called a type I point) or is < Con q (such a point will be called a type II point). 
Say that a point y is type III at time t if 2C n 9 < \ip y (t)\ 2 < an-2C n( 3+ «)/ 4 , 
and type IV if \ip y (t)\ 2 > an + 2C n( 3+( ')/ 4 . Note that if t g T, there cannot 
be any type III or type IV points. 

Let x be the mode at time 0. Since T, x is a type I point at time 0. 
Suppose x is a type III point at some time in A; let s be the first such time 
(there is a first time by the time-continuity of the wavefunction). Let s' be 
the last time before s such that x was type I at time s' . Then by continuity, 

an — 2Con^ 3+9 ^ 4 < \ip x (u)\ 2 < an — 

C on (3+?)/4 for all u e ^ }S y 

Thus, the interval (s',s) is a subset of T. Since |T| < e - n V4Co 5 j m plies 
that s-s' < e - n9 /4Co. Moreover, |^(s)| 2 < an-2C n( 3+( ?)/ 4 and \tp x (s')\ 2 > 
an — Con^ 3+<? ^ / ' 4 . Thus, there exists u G (s', s) such that 

C on (3+ 9 )/4 



(42) |l^WI 2 



- e -n 9 /4C* 



On the other hand by ([!]), for any y £ V and at any i, 
= 2Re(^y(t)^ y (t) 



V ryv 'dt 



Note that N(ip(t)) is the same for all t and is bounded by n c at i = (and 
hence for all t) by the assumption that g" T. Since \ip y (t)\ 2 < N(ip(t)) for 
any y and t and |2Re(it)| < 2\u\, this implies that for any y and i, 



< 



2|^(t)| 



/i 2 



^d^(t)i + i^(t)i) + i^(t)r 



< 4n 2p L>iV(^(t)) + N(ip(t)) 2 < n c . 

But this is in contradiction with ( |42[ ) if n is sufficiently large. The contra- 
diction is to the assumption that x is type III at some time in A. Thus, x 
cannot be type III at any time in A. 
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Similar arguments show that if \T\ < e ~ nq / AC o anc [ o T, and x is the 
mode at time 0, then x cannot be type IV any time in A, and for any y ^ x, 
y cannot be a type III point at any time in A. 

These three deductions combined with continuity of the flow and the 
conservation of mass and energy establish that x must be the mode at all 
times in A and moreover ([7]) must hold for ip(t) for all t £ A (with a C 
different than Co), thus completing the proof of the theorem. 



8. Proof of Theorem 2.5 (Limiting Distribution) 



Fix k distinct elements x\, . . . , Xk € V. Define an undirected graph struc- 
ture on £ as follows: for any two distinct a, r 6 S, say that (<r, r) is an 
edge if 

{ax%, ax k } n {txi, . . . , Tx k } ^ 0. 

Since £ is a group and Assumption 2 of translat ability holds, therefore 
for each a and each 1 < i < k, there are at most k automorphisms in S 
that take axi into the set {x\, . . . , Xk}- In other words, there are at most 
k possible r E S such that r _1 crxj € . . . , This shows that the 
maximum degree of the graph on S is bounded by k 2 . Consequently, there 
is an independent subset £' of S of size > n/k 2 . 

Let c^) = (</>x):rey be a vector of i.i.d. standard complex Gaussian random 
variables. Fix a Borel set t/CC*, and let 



W) := 



)eu}\ 



Since £' is a set of automorphisms and the Hamiltonian H is invariant under 
graph automorphisms, 

EQ(0)=P((</ w ..,^ fc )ElO- 

By the design of Q(0) is an average of independent random variables 
taking value in {0, 1}; consequently, by Hoeffding's inequality [20], for any 

e > 0, 

P(|Q(0) -EQ(0)| > e) < 2e" |E ' |e2/2 < 2e" ne2/2fc2 . 

8.1. Subcritical Limiting Distribution. If f3B 2 < 8 C , then as in the 
proof of Theorem |2.3[ for any e > 0, 

P(|Q(B-VV)-EQ(0)|>e) 

< e ~n"/c + e o^/ 5+Cn2 P +Cn29p(|Q( ^ _ EQW| > £) 

e -n q /C i 2 e Cn 4 / 5 +Cn 2 P+Cn 2 9-ne 2 /2fc 2 

A simple computation using the above bound and the identity 

/>oo 

E(X) = / F(X > t)dt 
Jo 
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that holds for any non-negative random variable X, shows that there is a 
constant C > depending only on the usual model parameters such that if 
n > C, then 

E\Q(B~ 1/2 iP) -EQ(<P)\ < kn- 1/c . 
By Jensen's inequality, this gives 

E\Q(B~ 1/2 ip) - EQ(<f>)\ > |EQ(B- 1 / 2 V) - EQ(0)| 

= IPCB-V 2 ^,...,^) e l/)-P((0 !Bl ,...,0 !Bfc )€^)|. 

This completes the proof in the case (3B 2 < 6 C . 

8.2. Supercritical Limiting Distribution. Next, suppose (3B 2 > 9 C . 
Take any q as in the first part of Theorem 2.2 and let Co denote the con- 
stant C. Let R denote the set of functions satisfying ([7]) with C = Co- For 
each x £ V, let R x be the subset of R consisting of functions with mode 
at x. Note that the R x s are disjoint and R = \J X R X . Let a = a((3,B) be 
defined as in @. Recall that 

F(J3, B) = Pa 2 + log((B - a)vre). 

Fix e > and let 

A := {/ : |Q((5 - a)' 1 ' 2 f) - EQ(0)| > e}. 
Then for n > Co, 

P(V> e A) < P(v» -R) + g A n i?) 

^e-^ + Z- 1 / e-W)l {S2(/)<Bn} d/ 
•/Ani? 

J Ann 

< e -n*/C + _ a)7re) -n e Cn^ + Cn^+Cn^ £ ^ r ^ 

Fix x £ V. For each /, let /' be the function such that = f y for all y ^ x 
and = n~ 1 / 2 f x . Note that the map / i— > f is linear with determinant n _1 
(and not n -1 / 2 , since we are contracting both the real and imaginary parts 
of/,). 

There are at most k many a G S such that a 1 x <E {x%, . . . , rr^}. Therefore 
for any /, 

|g((J3 - a)" 1 / 2 /) - Q((5 - a)-VV)| < ^ < ^- 

Consequently, ii f £ A then /' G A' where 

A' := {/ : \Q({B - a)^ 2 f) - EQ(<f>)\ > e - k 3 /n}. 
Again if / G R x , then /' G R 1 where 

R> ■= {/ : S 2 (f) < (B — a)n + Cn^ +q ^ 4 }. 



Combining the above observations gives 
Yol(AnR x ) = J l {feA nR x }df 

- n j l{f'eA'nR'}df' 

< ne Cn( ' i+q)/i+n j e-WViB-^tf 

= n{(B - a)7re) n e Cn(3+ " )/4 P(|Q(0) - EQ(<f>)\ > e - k 3 /n). 
The proof is now completed as before. 
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